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VELOCITY DISTRIBUTION IN INCOMPRESSIBLE POTENTIAL FLOW 
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SUMMARY 

A step-by-step numerical procedure based on a conformal -mapping 
theory is presented for the design of a cascade of airfoils with a 
prescribed dimensionless -velocity distribution in incompressible 
potential flow. The numerical procedure includes a set of tables to 
serve as a guide in computation. A numerical example is also 
presented in order to illustrate specific solutions to some of the 
more difficult- problems arising in the design of a cascade by this 
method. 

The conjugate harmonic function was evaluated by use of the 
method outlined by Naiman. Two separate calculations were made: one 
used 360 equally spaced values of the circle angle 0, and the 
other - used 90 equally spaced values of 0. The 90-point method did 
not insure sufficient accuracy for the example used. The accuracy, 
however, depends on the design conditions prescribed and, consequently, 
it cannot be inferred that the 90 -point method will give inaccurate 
results in all cases. • 


INTRODUCTION 

Several exact solutions to the problem of obtaining an airfoil 
cascade in incompressible potential flow with a prescribed velocity 
distribution on an airfoil In the cascade have been obtained by the 
method of conformal transformation (references 1 to 3) . The numerical 
computations involved in conformal -mapping methods are, however, 
us uall y quite long and complicated. A detailed numerical procedure 
that is intended to mahe the actual application of a conformal- 
mapping method to cascade design more efficient in regard to time 
and personnel was therefore developed at the NACA Lewis laboratory 
and is presented herein. 

The procedure employed requires the specification of the 
dimensionless-velocity distribution on an airfoil -in the cascade 
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an d the values of the upstream and downstream velocity vectors of 
the flow through the cascade. The velocity distribution specified 
for an airfoil may not he theoretically attainable; that is, the 
airfoil profiles calculated by use of this velocity distribution 
would not have closed contours. It is possible, however, to 
modify the velocity distribution in such a manner that closed pro- 
files will result. 

A numerical e xam ple is included with the procedure in order 
to illu strate specific solutions to some of the more difficult 
problems arising in the design of a cascade by this method. A 
brief s ummar y of the theory underlying the numerical procedure is 
presented to acquaint the reader with the fundamental principles 
employed in the solution of the problem. More detailed discussions 
of the theory are presented in references 1 to 3.- 


SYMBOLS 

The following symbols are used in this report: 

C oonstant 

d dimensionless spacing of airfoils in oaBoade (ratio of 

actual spacing to prescribed suction-surface arc length) 

f(0) function added to p(0) to satisfy conditions on h(z) 

h(z) analytic function of z 

k constant locating complex sources in z -plane . 

N constant determining amplitude of sine curve 

p(x,y) harmonic function 

p f (©) p(0) + f(0) corrected harmonic function 

Q(z) function of z defined by equation (11 ) 

q(x,y) conjugate harmonic function of p(x,y) 

Ee "real part of" 

dimensionless airfoil Eire length (ratio of actual Eire length 
to suction-surface arc length) 


S 
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S' 

^l' ®T,2 

u(S) 

u* 

Y 

V 
V f 

T C ( 0 ) 

V=iCf4-l\J/ 

z=x+iy, 

a 

P(0) 

r 

7 

A 

A k 

e 

lUS+iq, 

0 

% 


corrected dimensionless airfoil arc length (ratio of 
corrected arc length to prescribed suction-surface 
arc length) 

value of S corresponding to trailing edge of airfoil 
when approached from pressure and suction surfaces, 
respectively 

airfoil velocity function (equation (4)) 
corrected value of u 

magnitude of dimensionless velooity vector (ratio of 
actual velocity' to maximum airfoil velocity) 

magnitude of prescribed dimensionless airfoil velocity 
(ratio of actual velocity to maximum airfoil velocity 

corrected value of v 

velocity on unit circle 

complex potential function 

complex variable denoting circle plane 

constant determining phaseangle of sine curve 

function of 0 defined by equation (25) 

circulation (positive counterclockwise) 

constant deteuning wave length of sine curve 

velocity-potential range on airfoil 

velooity -potential range on unit circle 

included trailing -edge angle of airfoil 

complex variable denoting cascade plane (equations (23) 
and (24)) 

circle angle (z -plane) 
particular value of 0 
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<p velocity potential 

X angle of velocity vector (positive clookvd.se) 

A^ auxiliary constant defined "by equation (28) 

^ stream function 


Subscripts : 


a 

c 

1 

m 

U 

T 

1 

2 

3 

4 

I 


oascade plane 
circle plane 

indicates particular values of N, 7, and a 

mean flow conditions 

leading edge (nose) 

trailing edge (tail) 

upstream 

downstream 

various constants 


A function of (z = x + iy) or of (x,y) expressed as a funotion 

i Q 

of 9 Implies that the function is evaluated for z ■= e . 




REVIEW OF BASIC THEORY 

The solution to the problem of designing a cascade by the 
method to be outlined is baaed on conformal -mapping techniques 
similar to those employed in references 1 to 3. The desired cascade 
of airfoils is determined by forming an analytic funotion of a 
complex variable, which effects a mapping of the incompressible 
potential flow about a unit circle into the inco m pressible potential 
flow about the cascade. The mapping function is such that the uni t 
circle is thereby transformed into the airfoils of the cascade. 
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The numerical solution involves two problems: (l) the 
determination of the potential flow about the unit circle; enfl 
(2) the determination of the mapping function and. the airfoil 
coordinates . Before these problems can be solved, however, the 
following design information must be specified: 

1. The magnitude of the dimensionless velocity v along an 
airfoil in the cascade, as a function of the dimensionless 
airfoil arc length S 

The dimensionless arc length S will be considered negative from 
%,1 ' fco zero , and positive from zero to Sp g, vhere Sp and 

represent the trailing -edge stagnation point on an airfoil 

when approached from the pressure and suction surfaces, respectively 
The leading-edge stagnation points correspond to S = 0 (fig. 1). 

2 . The upstream and downstream dimensionless velocity vectors 
of the flow through the cascade 

These vectors will be designated Y^ e” i ^l “ s-nfl Y g e”^^2 ” n 
respectively, where Y^ and. Y 2 are the magnitudes of these vel- 
ocities and X^ and Xg are the directions (positive clockwise) of 

the velocities as measured from a line perpendicular to the cascade 
axis (fig. 2). The four quantities Yj, Yg, X-l, and Xg cannot 

be specified independently, however, as they must satisfy the 
‘continuity equation 


Y^ cos X^ » Yg cob X g 


Determination of Complex Potential Function in Circle Plane 

In this discussion of the solution of eaoh of the previously 
mentioned problems, the complex plane containing the cascade is 
designated the £= ( £ + ir|) plane or cascade plane; whereas the 
complex plane containing the unit oirole is designated the 
z =* x + iy plane or oirole plane. 

The complex potential function W (z) for the flow in the 

0 ±k 

circle plane due to complex sources located at z = ±e ( (ref- 

erence 4, notation modified) is given as 
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It 1 b evident that equation (l) Trill "be specifically determined 
when the constants k, V m , A m , d, and P are known. The con- 
stants 'V m ,A m , d, and P can he determined directly from the pre- 
scribed velocity distribution for the airfoils in the cascade and 
the upstream and downstream velocities of the flow through the 
cascade . Specifically, 


A 

m 


tan 


-1 / tan + tan A^ 

v 2 ; 


( 2 ) 


where 



where 


Tj cos A^ 


m cos A. 


m 


•?T,2 


u(s) dS 


u(S) = - v(S) for S^ = S ^ 0 
u(S) = v(S) for 0 < S = Srj^ 2 


(3) ' 

(4) 
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(v(S) is the prescribed -velocity as a function of S 


4 = 


T 

(v-j_ sin A 1 - Yg sin Xg) 


(5) 


The constant k is determined hy the condition that the range 
of the velocity potential (the difference between the maximum and 
minimum values of the velocity potential) on an airfoil must equal 
the range of the velocity potential on the unit circle. 


The velocity potential on an airfoil cp (s) can he determined 
from the equation a 


<t> a (s) 


S 

n 


u(S) dS 


T,1 


( 6 ) 


The velocity-potential range A on an airfoil is then 


A = ' c ^a (0) ( 7 ) 

The velocity potential on the unit circle z = e^® is given 



where the constant C 3 in equation (l) has been so chosen that 
<P C (0) is zero for the value of 0 corresponding to the trailing - 

edge stagnation point on the circle. Then the velocity -potential 
range on the circle A k for a given value of k is 


Ajj. = cp o ( 0ip + 2jt) - cp c (0jj) 


(9) 


where 0^ and 0jj are the values of 0 corresponding to the 

trailing- and leading-edge stagnation points on the circle, 
respectively . 

It is possible to find k (and also 0rp and 0^) by a process 
of trail and error so that A = A fc . (See section NUMERICAL PROCEDURE . ) 



Determination of Mapping Function and. Airfoil Coordinates 


O) 


!Elie general form of 
"between the £ -plane and 


a mapping function that effects the desired transformation 
the 2 -plane is given by 


d£ ae" k 
dz it 




e 2 V 



exp 



e k + 



( 10 ) 


where 

c specified included trailing -edge angle of airfoils in oasoade (for example, C = 0 for 
ousped tailB, c = it for rounded tails) 


Q(z) 


z 2 ( e 2k - e -2k) 
(e 2k - z 2 )(z 2 - e -2 ^) 


(ID 


and 


h(z) » p(x,y) + iq(x,y) 

is a complex function, analytic and single -valued for |z| > 1, and is such that 


lim zh(z) a 0 (12) 

z_yo 

It is seen from equation (10) that the mapping funotion is specifically determined if the 
analytic function h(z) = p(x,y) + iq(x,y) is known. Actually, the process of obtaining 
the coordinates of an airfoil in the oasoade requires only that d£/dz be known for values 
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of z on the unit oirole. It is therefore sufficient to evaluate h(e^®) m p(0) + iq.{0) 
in order to accomplish this result. 

Evaluation of p(8) and q.(fl) , - A condition that holds betveen the cascade plane and 
the circle plane is 


q> a (t) + (£) -w c («) 


vhere T a and '1^ are the velocity potential and the stream function, respectively, for 
the flow in the casoade plane. From this relation it is possible to deduo a that 


vhere 



v 0 (e) » Ka 
* % 


dV 0 (z) 

dz 
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(13) 


(14) 

(15) 


is the velocity on the unit circle and u(9) = u(s) as a function of 0. 
From equations (10) and (14) is obtained 


V q( 9 ) 

u(^) 


MSI 


2 sin 


( 9 ~ ^ 
\ 2 / 


- 1 - 1 


eip- 


Q 

V 


- h -I— - — 

* i2 2rt 


log^2 (coah 2k-*coa 2^,) - 


k ( cos 0, 

k - e cos 9 tanh 


9 ? . 

cosh k / ysinh k, 


(16) 


<o 
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Equation (16) can now "be solved for p(0), 

M 


p(e) = Q(e) 


f * 

l log e Wey ~^ej 


2 sin (LZ JH I 


± - 1 
It 


k + (i - ~ 2a) log e \/ 2 ^ cosh 2k “ cos 20 T ) “ k - 


k a . ,-l /cos 0m \ fc . . . _1 /sin ftpS 

e cos 9 tank I . — — ) - e sin 0 tan A 

\cosh ky \sinh k> 


(17) 


•where 


jc < 
" 2 “ 



< £ 
= 2 


All quantities in equation (17) have keen previously evaluated 
except u(0). The procedure for evaluating u(0) is as followBi 
Erctn equations (8) and (13) the relation 

«P a (S) = q> c ( 0 ) (18) 


'is obtained for corresponding values of S and 0. By matohing 
these velocity potentials, S can he obtained as a function of 0 
and. consequently, u(s) oan be obtained as a function of 0, 
u(0). Hence, p(0) can be evaluated from equation (17). 


The necessary and sufficient condition that the condition 
Imposed on the function h(z) by equation (12) 1 b satisfied is 


0 T + 2jt 


0 T + 2« 


p(0) d0 ** 

a p 

p(0) sin 0 d0 = 

0m 

0m Je 


0rp + 


T 


2ir 

p(0) cos 0 d0 


0 


Usually the values of these integrals will not be zero and, 
consequently, a correction term f (0) must be added to p(0) such 
that if 
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p'(e) B p(0) + f(8) 


Then 


+ 2it 


p' (0)40 


0 T + 2it 


0q> + 2jt 


J0n 


p' (0) sin 8 40 


T 


p’(0) cos 0 40 = 0 


(19) 


If it is assume! that the values of the integrals involving 
p(0) are not zero, "but instea! sere 


0qi + 2it 


p(0)40 = C 4 


0fp + 2rt 


p(0) sin 040 o C c 


T 


0^i + 2it 


p(0) cos 040 = Cg 


T 


it follows that f (0) must "be a function such that 

O 

2 * 

f (0)40 = -C 4 


0 T + 2* 


0TJ1 

0tj + 2it 

f (0) sin 040 = -C 5 


* 


0 T + at 

^ f(0) cos 040 = -C, 


T 
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The proper choice of f (0) fulfilling these requirements is 
somewhat arbitrary. The addition of f(0) to p(0) will, however, 
produce a change in the prescribed airfoil velocity and the pre- 
scribed arc length on the airfoil for 0 T = <0 £ 0 T +2 rt, given by 


and 


v'(0) 

v'(0) 


- u(0)e _f ( 0) / Q(0 ) for 0|£ = 0 = 0jj 
u(0)e" f ' for 0 N = 0 = 0 T + 2rt| 


( 20 ) 


S’(0) 


0 

r 


e 


f(s)/<i(0) 


% 


d0 

d0 


( 21 ) 


These .equations are the parametric equations for v’(S') 
where v'(S') represents the corrected airfoil velocity as a 
function of corrected arc length. It is possible in many instances, 
however, to prescribe f (0) in such a manner that it is zero for 
values of 0 corresponding to certain regions of the airfoil, such 
as the suction surface. Over these regions, the magnitude of the 
prescribed velocity as a function of 0 will remain unaltered. 


The function q(0) can be evaluated for various values of 0 
from the equation 


Jtji + 2rt 


q(0 i ) = 


p(0) cot (- 


0 - 0 , 


-) d0 


( 22 ) 


where 0^ is a particular value of , 0. The. function p(0) in 

this equation is replaced hy p'(0) if p(0) has been corrected. 

\ 

Determination of airfoil coordinates . - The coordinates (£ , T|) 
of an- airfoil in the cascade can be obtained in the following manner j 
T he complex variable £ = | + iq can be expressed as 



1242 



1242 


MCA TN 2101 


13 


The real and imaginary parts of &z for z = e ±w now give 

Jdz 

the coordinates £ and tj, respectively, of the airfoil. Specifi- 
cally, for p(0) uncorrected. 


.10 


t = 


dt 

dz 


cos P (9) dS = 


Y (0 ) 

-41 cos p * 9 


(23) 


T) = 


dt 

dz 


sin P (0) d0 = 


y (0) 

_2__ sin p (9) d9 

u(0) 


(24) 


where 


dt> 


P(9) =|. + 0 + arg (g) i0 

Z —0 


e 

2 


it _ 0 ^ l + .«-i r fcan e a 


=. aM + (i -£)<*-£ + i tan- 

0(01 JC )2 2 2 


vtanh k/ 


cos 9 tan 


. / sin 9rp\ 

-J. _ sin 0 tarih' J - i r— r- 

V siuh V \ COBh 


^ / cos 0fj^ 


+0+ | . (25) 


_n ( tan 0m\ 

and where tan \ — is in the same quadrant and has the same 

\tanh k/ 

sign as 9^. 

If the function p(9) has "been corrected, the following 
equations are used in place of equations (23) and (24) : 


v^ 00B p (9) as 


(26) 


ij = 


Za^sin p (8) ae 


(27) 
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Ee strict ions on Airfoil Velocity far Specified Trailing-edge Angles 

It can "be shown from equation (17). that the function u(0) 
most possess the following properties for cusped-, rounded-, or 
pointed-trailing -edge airfoils: 

1. Cusped trailing edge: - u(© T ) = +u(0 T + 2it) = C, where 

C is a positive constant 

2. Bounded trailing edge: u(© T ) = u(0 T + 2 it) = 0, where the 

zero is of order of one 

3. Pointed trailing edge (trailing-edge angle = e ): u(0^) «= 

u(0 T + 2it) = 0, where the zero is of order of e/rt. 

As there is no convenient method of determining from u(S) the 
behavior of u(0) at ftp and 0^ + 2rt before extensive calculations 

have been performed, it may be necessary to modify u(0) to fulf i l l 
the preceding conditions . This modification will result in a slight 
change in u(S) at S^, ^ and S^, 2 . An approximation to the pre- 
ceding conditions on u(0)- can be made, however, by applying the 
restrictions to u(S) at ^ 1 and 

NIMEHICA1 PBOCEDUEE 

The numerical procedure that follows will liet the computational 
steps leading to the evaluation of the coordinates of an airfoil in 
the cascade . These steps have been included in a set of tables to 
serve as a guide in performing numerical computations. The two 
problems previously mentioned - namely, the determination of the 
complex potential function in the circle plane, and the determi- 
nation of the mapping function and the airfoil coordinates - will 
be treated separately. 


Determination of Complex Potential Function in Circle Plane 

The complex potential function W(z) (equation (1)) can be 
determined if the constants ?v m , V m , T, d, and k are known. 
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Evaluation of Xjjj, "^in> F* and d. - The constants 

V , T, and d are evaluated from equations (2), (3), (4), and 
(5), respectively. 

Ev a lu a tion of k . - The constant k is evaluated "by a process 
of trial and error that simultaneously determines 0^ and 0^. 

The procedure required to evaluate these constants is: 

1. Determine the velocity-potential range A on an airfoil in the 
cascade from equation (7). 

2. Assume a value of k (probably between 0 and 1) and determine 
the auxiliary constant 


where 


= tan 


-1 


(tan ^ tanh k) 



3. Determine 0^ and 0^ from the following equations, 
(reference 4) : 


0 N + \ = sin -1 


r -r 

vW 


cosh k sin 
8inX m 



(28) 


(29) 


ftp = - 0jj - jt - 2X^ 


(30) 


4. Calculate the velocity potential range on the unit circle for 
the assumed value of k from 



16 


KA.CA TN 2101 


A k =cp c (0 T + at) - cp c (e N ) 


(9) 




tan 


_3_ /sin 0ip N 


. sink k ; 


- tan' 


-1 f Bin V 

, sink k> 


cos A 


m 


/cos Sm\ , , /cos 

tank ± I , — _ ] - tank x ‘ ’ 


.cosh k. 


.cosh ky 


where 


2T m cL 


tan 


f tan ftp 
\tank k 


+ 2it - tan' 


/ tan ft 


vtank 



- — <tan 
2 


- 1 ( sin g T \ <1L 
vsinh k/ 2 


2 \sinh k/ . 2 


-i 

1 \tanh t) 


and tan 
and have the same sign as 


* + -1 / tan0 n\ 
**■ taa vt^kT; 

and Sjj, 


ft 


are in the same quadrant 
respectively. 


5. Bepeat steps (l) to (4), assuming several nev values for k 
until, ky interpolation or plotting against k, a value of 

The correct values of k, 9^, 


k can he found such that 


A k = A. 


a nd ftjj are those that produce this last equality. (The evaluation 
of k is given in table I.) 


Determination of Mapping Function and Airfoil Coordinates 

Tke mapping function -will he determined on the. unit circle 
if the functions p(0) and q (ft) are known. The airfoil coordinates 
can then he obtained from the mapping function. 
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Determination of -p(0). - 

1. Evaluate the velocity on the unit circle for (% = 0 % 0™ + 2it) 

, ri 

d¥.(z) 


v„(e) = Ee 


dz 


z=e 


±e 


(15) 


_ Y^d sin X^ Q(0) _sin (0 + X - sin (0jj + X^)J 
- — ; 

where Q (0) is given "by 


sin X^ cosh k 


g/ 0 \ _ sinh 2k sinh k cosh k 

cosh 2k - cos 20 sihh^ k sin^ 0 

(The evaluation of v c (0) is given in columns . 1 to 9 in table II) 

2. Evaluate u(0) in the following manner: 

Plot q? a (S) (equation (6)) and evaluate cp c (0) from 
0^i = 6 ^ 0q, + 2 it, where 

<P C (0) = Ee [^(e 10 )] 


( 8 ) 


= ^2— ( | sin X„ tan -1 ( - cos X TO tanh -1 ( ??£, 


fr 


* 


m 


,-1 / sin 0 ^ 
\^sinh k^ 


'm 


.cosh k/ 


_ .t.an-1 / tan 0 + c 

2Y m d \ tanh k ft 7 


where Cy is a constant such that cp Q ( 0^ ) = 0 and 


- *<tan \ <* 
2 l^sinh k ) 2 


The term tan”-'- ( ^ . an ... - ... ) is in the same quadrant and has the same sign 
\tanh k J 

as 0. (The eval u ation of <P Q (0) is given in columns 10 to 24 in ■ 
table H.) 
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Beginning at the point ^ on the airfoil velocity -potential 
curve and at the point 0^ on the circle velocity -potential curve, 
progressively match the values of cp a (S) and. cp Q (0) from 
and 0 T to Sj 2 ani %> + 2jt * thereby establishing a one-to- 
one correspondence between the values of 0 and S, which satisfies 

<P a (s) »cp o (e) (is) 

As S(0) is known, u(S) is obtained as a function of 0, u(0) . 

3. The function p(0) is now determined from equation (17) (The 
evaluation of p(0) is given in columns 25 to 38 of table II.) 

Correction of p(0) to satisfy conditions on h(z) . - Compute 

rrals 

d0 


sin 0 d0 


cos 0 d0 

as indicated in columns 39 to 43 in table H . 

If these Integrals are not zero, form the function p’(0) = p(®) + f(®) 
where p'(0) and f(0) fulfill the conditions of equation (19) 
(columns 44 to 50 of table H) . One expression for f(0) is given by 

f (0) = ET-l sin 7 ± (0 + c^) + Ng sin (0 + ag) + N 5 sin 7 3 (® + &$) 

(31) 


(by numerical Integration) the irrfceg 


0m + 2jt 

n 

P(0) 

0m 


n 


p( e ) 


0 T + 2jt 


P(9) 
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for 


■where 3 = 1,2,3 
to J; f(0) = 0 

The constants N^, 7y and (j = 1,2,3) -are chosen so 

as to fulfill the integral requirements and, if possible, restrict 
the velocity change to a minimum. It seems advisable in this case 
to assume various values for 7. and a. and to evaluate N, from 
the integral requirements. J J 

Determination of q(0). - The function q(0) is evaluated frcm 
equation (22) either by some form of numerical integration or by the 
method of CTaiman (reference 6) . 

Determination of airfoil coordinates . - 

1. Evaluate 0(0) from equation (25). (The evaluation of 
0(0) is given in. columns 51 to 64 in table II.) 

2. Evaluate the airfoil coordinates from equations (23) and (24) 
if p(0) has not been corrected, or from equations (26) and (27) if 
p(0) has been corrected. (The airfoil coordinates are evaluated in 
columns 65 to 71 in table H.) 

APPLICATION OF NUMERICAL PROCEDURE 

Several important problems arise in the application of the 
numerical procedure to the design of a particular cascade. These 
problems may be listed as follows : 

1. The number of values of 0 to be chosen for computations 

2. The method of selecting functions to add to p(0) in order 

to satisfy conditions on h(z) 

3. The method of evaluating q(0) from equation (22) 

4. The type of numerical integration to be used in evaluating 

the various integrals that appear in the procedure ' 


(- v -iL)Sa S- aj 


in the term that has the subscript corresponding 
for all other values of 0. 
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The solutions to these problems of the numerical techniques 
to he employed to give a certain aocuracy in computations will 
depend largely on the particular design being considered. At present^ 
no general solutions can be given. 

In order to illustrate particular solutions to these problems, 
however, and to clarify the numerical procedure, a cascade of air- 
foils with cusped tails was computed. It was decided in advance to 
compute an airfoil in the cascade in two ways. One method was to 
use 360 equally spaced values of 0 and the other was to use 90 
equally spaced values of 0. The choice of the number of values of 
9 was made in order to investigate the accuracy of using the method 
of Naiman (reference 6) in computing the values of q (0) for the 
problem chosen and to investigate further the influence of this 
accuracy on the computed airfoil coordinates. The values of 0 had 
to be equally spaced to use Naiman's method for computing q(0). 

Naiman' s method had particularly been chosen for the evaluation of 
q(0) because of its applicability to calculation by punch-card 
machines. In all calculations where numerical integration was required, 
Simpson's l/3, or parabolic method, was used because it was felt 
that this method would give sufficient accuracy for the values of 0 
chosen. In all numerical work, five significant figures were used, as 
it was believed this accuracy was sufficient for the example. 


Design Specifications 

The velocity on an airfoil in the cascade was expressed non- 
dimensionally as the ratio of the actual airfoil velocity to the 
maximum airfoil velocity. The airfoil arc length was -expressed non-, 
dimensio n all y as the ratio of the actual airfoil arc length to the 
total airfoil suction-surface arc length. The velocity distribution 
is shown in figure 3. 

The magnitudes of the inlet and exit velocity vectors were chosen 
as, respectively. 


Y 1 = 0.68229 
Y 2 = 0.60000 

The angles of the inlet and exit velocity vectors were chosen as, 
respectively 


A 1 » 30.0000° 

A 2 = 10.0000° 
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360 -Point Airfoil Design 

The first step in the design of the 360-point airfoil in the 
cascade was the calculation of the velocity potential distribution 
by numerical integration of u(S) . Simpson^s parabolic rule was 
used taking values of S in intervals of 0.025. The velocity 
potential distributions cp a (s) and u(S) are shown plotted jointly 

in figure 4. From the velocity potential, T was determined to be 
0.28201 and cp^^, - 0.51556. The constants X , V m , and d were 

calculated from equations (2), (3), and (5), respectively. The values 
obtained were as follows: X m = 20.6484°, V m = 0.63144, d = 1.19012. 

The values of k, X k , %, and calculated by the trial- 

and -error process indicated in the numerical procedure, were found to 
be k = 0.48337, X k = 9.6614°, 0^ = - 15.6038°, and = - 183.7190°, 

respectively. 

The velocity distribution v Q (0) and the velocity potential dis- 
tribution cp c (0) on the unit circle were calculated from equations 

(15) and (8), respectively. These distributions are shown plotted in 
figures 5 and 6, respectively. 

The harmonic function p(0) was' calculated from equation (17) 
and is shown plotted in figure 7. 

In order to determine u(0), the plot of cp & (s) was superimposed 
upon the plot of u(S) (fig. 4). Each value of cp c (0) (fig. 6) was 
then located on the <P a (S) curve and the u(S) value was read for 

the corresponding value of S. In this manner, the value of u(S) 
for a given 0, u(0), could be directly determined. 

In order to determine whether or not p(0) fulfilled the con- . 
dltions on h(z), the integrals 


+ 2it 

f p(0) d0 
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Sfj + 2jt 

p(0) sin 0 d 9 


9 ^ + 2rt 
n 


p(0) cos 9 d 9 


1 


were numerically integrated and the following values were obtained: 

2 it 

p(0) d0 = 0.38236 


0rji + 2lt 


+ 2* 


*p (0) sin 0 d0 « - 0.09526 


T 


0 T + 2i( 


p(0) cos 0d 0 = 0.26919 


0„ 


In order to correct for airfoil closure, a function f(0) was 
formed such that 


0qi + 2jt 


0rji + 2rt 


0qi + 2jt 


J 9, 


p* (0) d0 = 


p' (0) sin 0 d0 = 


T 


p‘(0) cos 0 d0 <= 0 


( 19 ) 
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where 

I 

p*(S) = p(0) + f(0) 

Initially, f (0) was assumed to "be a function of the type 
given "by equation (31) . The wave lengths and zeros of the individual 
sine terms were picked somewhat arbitrarily except for the restriction 
that corrections on p(0) were to occur principally on the pressure 
surface. The amplitudes of the sine terms were then determined by 
the s imu ltaneous solution of the equations 

®rp + 2lt 0rp + 2K 

p(0) d0 = - 0.38236 


2k 

p(0) oos 0d© = - 0.26919 ' 


2k 

p(0) sin 0d0 = 0.09526 


Several other choices for the wave lengths and zeros of the sine 
terms were then made and the amplitudes for the terms were computed. 
The wave lengths and zeros that gave relatively small amplitudes far 
the terms were prescribed for the f(0) used for the correction. 

The function f(0) was defined in the interval 0^ 4 0 4 0 T + 2 k 

as follows : 


f(0) d0 


T 


T 


0<ji + 2 k 


+ 


f(0) cos 0d0 = - 


+ 2k 


0J + 


f(0) sin 0d0 = - 


T 


T 


for 


f(0) = 0 

- 183.719° = 0 < - 176.281° 
f (0) = - 0.21488 sin 7.5 (0 + 176.281°) 
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for 

- 176.281° = 9 < - 157.5° 

f(0) = - 0.21488 sin 7.5 (9 + 176.281°) + 0.008703 sin 4/3 (0 + 157.5°) 
for 

- 157.5° = 0 < - 152.281° 

f (0) = 0.008703 sin 4/3 (0 + 157.5°) 

for 

- 152.281° = 0 < - 27.719 

f(0) - 0.008703 sin 4/3 (0 + 157.5°) - 1.26792 sin 7.5 (0 + 27.719°) 
for 

- 27.719° = 0 < - 22.5° 

f(0) = - 1.26792 sin 7.5 (0 + 27.719°) 


for 


22.5 = 0 < - 3.719 

f (0) = 0 


for 


- 3.719 = 0 = 176.281° 

Graphically, f (0) plotted in figure 7 , consisted of the sum 
of three separate sine loops. The curve p’(0) = p(®) + f(®) is 
shown plotted in figure 8. 

The corrected velocity v’(0) and the corrected arc length 
S'(0) were computed from equations (20) and (21), respectively. The 
corrected-velooity distribution v’(S') is shown plotted in figure 9. 
This velocity distribution is the actual distribution for the computed 
airfoil shape. The peak on the suction surface in the nose region is 
undesirable, but steps were not taken to correct it in the present 
example . 
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The values of p* (0) at the 1° intervals were used to calculate 
the conjugate harmonic function q(0) (fig. 10) by use of the method 
outlined by Caiman (reference 6) . The actual computation of this 
function ■was performed on punch-card machines. ■ Interpolation was 
.necessary to find the values of q(0) as the preceding method gave 
values of q(0 - l/2°) . 

The values of q .(0) vere used to compute 6(0) from 
equation (25) and the coordinates of an airfoil in the cascade from 
equations (26) and (27) (numerical integration) . The resulting 
airfoil i-; shown in figure 11. 


90-Point Airfoil Design 

The computations for the 90-point design remained the same 
as in the previous case, in that the choBen 90 points were taken 
from, the 360 points except that the function p(0) had to be 
recorrected for closure. Eecorreetion was necessary because of 
the discrepancy in the numerical integration used, as the integration 
was performed over larger intervals » The corrections for closure 
were made within the limits chosen in the previous case. The result- 
ing function p'(0) is shown in figure 12. 

The conjugate harmonio function q(0) was computed from 
p'(0) using the 90-point method of ITaiman (reference 6) and is 
shown in figure 13. The airfoil coordinates were computed and the 
resulting- profile Is shown in figure 14. it is seen that this air- 
foil Is crossed near the trailing edge and remains open at the 
trailing edge. 


Comparison of 360- and 90-Point Airfoils 

It can he seen frcrn equations (25), (26), and (27) that a 
change In q(0) would introduoe a change in the values of the 
integrands in equations (26) and (27) and, consequently, would 
introduce a change in the final airfoil coordinates. Hence, the 
fact that the 90-point airfoil remained open although the conditions 
for closure were satisfied could possibly be attributed to the fact 
that the 90-point method of Hainan did not insure sufficient accu- 
racy in the values of q(0) for the present example. This pos- 
sibility was therefore investigated. 

In order to increase the accuracy of q(0), it was decided to 
use the 360-point method of Haiman. In order to obtain the 360 vali es 
of p’(0) necessary for the computation, parabolic interpolation was 
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used. This type of interpolation -was used "because the integration 
of p'(0) "by Simpson's 1/3 rule in satisfying closure requirements 
had already assumed a parabolic variation between the values of 
p'(6).' The interpolated values were then used to derive new values 
of q(0) at the chosen 90 points, and the airfoil coordinates were 
then computed o The new profile, shown in figure 15, has a closed 
contour. The slight variation between the shape of this airfoil 
and the original 360-point airfoil (fig. 11) can be attributed to 
the integration of equations (26) and (27) over larger intervals and 
the recorrection of p(0) for closure. 

An investigation was also made to determine, if possible, in 
•what manner and in what regions the change in q(0) influenced the 
Integrands of equations (26) and (27). In order to make this 
determination, the following calculations were made and plotted: 
the difference in the values of q(0), as obtained from the 
360-point and 90-point methods of Kalman (fig. 16); the differences 
in q(0) divided by Q(0) ; and the dif f erences ' in the integrands 
of equations (26) and (27); as obtained by using the values of q (0) 
derived from the two different methods. A plot of the differences 
in the integrand of equation (26) -was superimposed upon a plot of 
the difference in q(0)/Q(0) (fig. 17). A similar plot was made for 
equation (27) (fig. 18). 

The differences in q(0)/Q(0) show a definite correlation to 
the differences in the integrands in equations (26) and (27), as 
q(0) enters into the expression for 6(0) as q(0)/Q(0). It should 
be noted that the changes in- q(0)/Q(0) seem to have the least 
influence on the integrands in the vicinity of 0^ (- 15.6038°). 


■ . CONCLUSIONS 

It can be concluded from the comparison of the 90- and 360-point 
airfoils that the 90-point method of Naiman does not insure sufficient 
accuracy in values of the conjugate function q(0) for all cases, as 
shown by the numerical example presented. The accuracy of Naiman' s 
method will depend on the shape of the harmonic function p(0), which 
in turn depends on the design conditions prescribed. As a consequence, 
it is not to be inferred that Naiman 1 s method, applied to 90 points or 
less, will not give accurate results for most cases. 


Lewis Flight Propulsion Laboratory, 

National Advisory Committee for Aeronautics, 
Cleveland, Ohio, September 15, 1949.- 
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TABLE I - 00TLIKE FOR COMPUTATION OF k 



Definition 

Operation 

Remarks 

k . 

Assumed 


1 

sinh k 

sinh k 


2 

cosh k 

cosh k 


3 

tanh Tc 

tanh k 


4 

sinX m 

sinX m 


' 5 

COB X m 

cos X m 


6 

tan A ffl 

tan X m 


7 

tan Xt 

(6) x (3) 


8 

x k 

tan-1 (7) 

2 * Xk * 2 

9 

sin Xfc 

sin (8) 


10 

sin (eg + X k ) 

[ifer* < 2 > -(«]/(« 


11 

% + Xjj- 

sin-1 (io) 

2 + Xk»| 

12 

% 

(H) - (8) 


13 

0T 

- (12) - 2(8) - 180° 


14 

sin 0jj 

sin (12) 


15 

cos % 

cos (12) 


16 

tan ©jj 

tan (12) 


17 

sin ftp 

sin (13) 


18 

cos Qj 

cos (13) 


19 

tan 

tan (13) 


20 

sin ©j/sinh k 

(17)/(1) 


21 

sin ©ij/sinh k 

(14)/(1) 


22 

cos ©j/oosh k 

(!8)/(2) 


23 

oos ©jj/cosh k 

(15)/<2) 

■■■ 
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TABLE I - OUTLINE FOR COMPUTATION OF k - Concluded 



Definition 

Operation 

Remarks 

24 

tan 0j/tanh k 

(19)/(3) 


25 

tan ©if/ tank k 

(l6)/(3) 


26 

1 (sink Te) 

tan - ^ (20) 

Express in radians 
X < tan -l f Bln *1) < K 

■ 2 “ tan Isiartv = ? 

27 

^ M) 

tan - ^ (21) 

Express in radians 

- R . £ tan” 1 ^ 2 21 

2 Vsinh k/ 2 

28 

- /COS ©rp \ 

tafih ' 1 fearsJ 

tanh- 1 (22) 


29 

‘—1 &&) . 

tanh~l (23) 


30 

— &£) 

tan -1 (24) 

Express in radians; same 
quadrant and sign as 0j 

31 

/taa%\ 

^ tears; 

tan -1 ( 25 ) 

Express in radians; same 
quadrant and sign as % 

32 

r / sin %\ 

sin X m [tan-J- ^-fgrk/ 

-- 1 er?3 

(4) [( 26 ) - (273 


33 

T . 1 ( COB 

COB X m [tanh- 1 -^ J- 

t ( C0B %\1 

taDh VcostT/J ‘ 

( 5 ) [(28) - (29)] 


34 

r T -if tane r\. 

2V m d 'tanh k/ 

tan_1 (Sbte) + 2n ] 

20 & 30 ^" 
m 

(31) + 6.2832] 


35 


(32) - (33) + (34) 

' 

36 

A k 

•V M d 

-r * (3 5) 

This result must equal 
F-<Pnri n for correct value 
for k 
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TABLE H - OOTLXHE FOR COMPUTATION OF AIRFOIL COOBDIHAIES 



Definition 

Operation 

Remarks 

e 

Assign in e<jual intervals 

0 Q 1 m 0 ■ Qrji + 2?l. 

1 

26 

29 


2 

cos 26 

COB (1) 


3 

cosh 2k - cos 29 

cosh 2k - (2) 


4 

. . sinh 2k 

^ cosh 2k - cos 26 

sinh 2k 

~ XT) 


5 

e + x fc 

9 + Xfc 


6 

sin (6 + Xj;-) 

sin (5) 


7 

sin (6 + Xfc) - sin (% + X^) 

(6) - sin (Ojj + X k ) 


8 

T T m d sin X m 1 
it sin X k oOkh k 

[sin (6 + Xjj) - sin (6u + X k )] 

— * ( 7) 
n sin cosh k 


9 

v c (e) 

(8) x (4) 


10 

sin 9 

sin 6 


11 

cos 6 

cos 9 


12 

tan 6 

tan 6 


13 

sin e/sinh k 

(10)/sinh k 


14 

cos 6/cosh k 

(ll)/coBh k 


15 

tan 6/tanh k 

(12)/tanh k 


16 

tan_1 (sihrr) 

tan-1 ( 13 ) 

Express in radians 

- £ S tan-l i £ 

2 \sinh k / 2 

17 

—- 1 (Srs) 

tanh”! (14) 
1 


18 

--MSte) 

tan"l ( 15 ) 

Express in radians 

tan ' 1 (estt) 1s in 

the same quadrant and 
of the same sign as 6 
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TABLE H - OOTUUE FOR COMPUTATION OF AXRFOIL COOBBITTATKS - Continued ' 


Definition 


Operation 


Remarks 



C 1 - k - (33) - (34) 0^^ = logg /2( cosh 2k - cob 20^, 


p(g) ae 


&r+2rt 

r 

/ p(6) si 11 6 d8 


p(e) cos 6 ae 


45 p'(0) « B(e) + f(8) 


p'(0) sin 0 


p'(0) cos 0 


p'(0) sin 9 d0 


e = Jt for rounded tail 
e n 0 for ousped tail 


(2 ” arc) * ^ 


(32) + k + (36) 


(4) x (37) 


(38) x (10) 


(38) x (11) 


(38) d9 


(39) ae 


(4o) ae 


Assigned arbitrarily] See numerical Procedure 




(38) + (44) 


(45) x (10) 


(45) x (11) 


(45) ae 



Must te equal to zero 


Must Le equal to zero 


\NACA^ 
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TAELS H - OUTLINE FOE COMPUTATION OF AIRFOIL COORDINATES - Continued 










































Airfoil velocity, v(S) 
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Figure 3. - Prescribed velocity distribution on airfoil. 
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Airfoil velocity potential, 9 a (S) 
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Figure 5. - Velocity distribution on unit circle. 
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p(0), f(0) 



Figure 7. - Harmonic function p(e) and correction function f(9), computed for 360 points. 
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Figure 8. - Corrected harmonic function p'Cfl), computed for 360 points. 



Modi fled velocity, v' (S') 
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-183.719 -143.719 -103.719 -83.719 -23.719 16.281 36.281 96.281 136.281 176.201 


C t 1 * q n 1 a D Ann 

w ^ 1 v ^ ‘ y ' u 


Figure 10. - Conjugate harmonic function q(0), computed for 360 points. 





Corrected harmonic function, p’(6) 
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183.719 -143.719 -103.719 -63.719 -23.719 (6.281 56.281 96.281 136.281 176.281 

Circle angl e. 0. deg 


Figure 16. - Differences between 36Q-polnt q ( 9 ) and 90-polnt q(0) 




Aq ( 0 ) /Q ( 0 ) , Ai 



-183,719 -143.719 —103- 719 -63,719 -23,7!9 !6,28! 56,28! 96,28! (36=28! !76,28! 

. Ct rcl e angl e, 8, deg 

Figure 17.- Dl fferences In 360-point q(9) and 90-potnt q(9), divided by Q(9); 
differences In t Integrand from 360-potnt q(9] and 90-polnt q(0). 
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1 83. 719 -143.719 -103.719 - 63.719 -23.719 16.281 56.28 ! 96,281 136.281 176.281 

Ci rcl e angl e, 8, deg 

Figure 18, - Differences In 360-point q(9) and 90-polnt q ( 9 ) , divided by Qle); 
differences in n integrand from 360-point q ( © ) and 90-polnt q|Sj. 








